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ABSTRACT 


The dynamical, chemical, and thermal evolution of zero-metal gas clouds 
has been modeled to study conditions of star formation in the early universe. 
Numerical results are given for the collapse of spherical clouds of mass 10^ 
and 5 x 10 r M G . Cooling by H z lines and by photons emitted in H + e~-*H- + 
h^ maintains collapse until formation of an equilibrium protostellar core of 
mass 0.02 M e . The cooling by photons produced with H is essential for low- 
mass star formation. If the cloud is fragmented, the evolution of the pieces 
is similar to that of the parent cloud, but the equilibrium core has larger 
density and mass. 



I. BACKGROUND 


The objective of the work performed under this NASA/Ames-UCSC 
Cooperative Agreement has been to numerically model the collapse and cooling 
of gas clouds containing hydrogen and helium but no heavier elements. This 
composition is chosen to represent the early universe following Big Bang 
nucleosynthesis, for which theoretical chemistry calculations predict a 
composition of 75 percent hydrogen and 25 percent helium with only a small 
fraction (/vlO - **) of heavy elements, i.e. ''metals" (Wagoner 1973). The basic 
questions addressed by our work are whether stars can form as gravitational 
condensations within this gas and whether their masses and evolutionary 
histories differ from those of objects with the higher metal abundance 
characteristic of our galaxy today. A fundamental problem has been that 
stars completely devoid of the heavier elements are not seen: the smallest 
metal abundance observed in our Galaxy and its globular clusters is about 
10 -5 \ This has led to speculation that the first generation stars were all 
sufficiently massive (>1 M 0 ) to have completed their evolution by the 
present epoch, and therefore it is of interest to inquire whether there are 
theoretical arguments favoring the formation of massive stars. 

A lower limit to the mass of gravitationally unstable perturbations in 
a homogeneous, static medium is the Jeans mass My , which varies with density 
and temperature as Mj-^T^/n^ 1 - (Jeans 1928). A commonly employed model for the 
formation of galaxies and stars is the process of hierarchical fragmentation 
in which the background medium forms unstable condensations of progressively 
smaller size (see, e.g., Hoyle 1953). The gravitational energy of a collapsing 
condensation is converted ultimately to- thermal energy, and if there is no 
means for cooling the gas, the temperature will increase as n x/3 , and the 
Jeans mass will increase as n t/z . Thus for continued fragmentation the gas 
must radiate sufficient energy that the temperature increases more slowly than 
n ,(3 . This may occur through collisional excitation of atoms or molecules 
followed by photon emission, but the excitation energies of atomic hydrogen 
and helium are so large that these species cannot cool the gas until the 
temperature approaches 10 M K. If, however, a sufficient fraction of the 
hydrogen exists as H a , the gas can radiate at and therefore maintain a much 
lower temperature through the excitation of low-lying rotational and vibra- 
tional states. 

It is widely assumed that molecular hydrogen is formed on dust grains 
within the dark clouds of our Galaxy, but it is difficult to envision grains 
in a gas without heavy elements. Mechanisms for forming H z in the 
primordial gas must be quite different from those acting today. The direct 
radiative association 


H + H — > + hV (1) 

is too slow for significant H*. production within the present age of the 
universe (Gould and Salpeter 1963), but residual ionization from the Big 
Bang may have allowed production by ion-neutral reactions. Peebles and 
Dicke (1968), Takeda, Sato, and Matsuda (1969), and Hirasawa, Aizu, and 
Taketani (1969) showed that the reactions 

H + e” — => H~ + hv (2) 

H"+ H — > H z + e” (3) 

yield an fraction of 10 during collapse to density 10 cm' J . At higher 
densities Palla, Salpeter, and Stahler (1983) find that the three-body 



reactions 


H + H + H — =9 H z + H (4) 

and H + H + H 2 -=> H x +H z (5) 

eventually convert almost all of the hydrogen to molecular form. 

Prior to recombination, density fluctuations in the expanding universe 
are dissipated by radiation pressure," but following the decoupling of matter 
and radiation at T-^-4000 K, n^lO* cm" 3 , bound clouds can form and contract 
to high density (Peebles and Dicke 1968). Hutchins (1976) has modeled density 
perturbations which expand with but more slowly than the background universe 
for a time t w and then collapse. Both the perturbation and the background 
are assumed to evolve as Friedman models, characterized by a parameter 
which determines whether the model is open, flat, or closed. Standard cosmo- 
logical equations are solved for the background density and temperature and 
the density of the perturbation as functions of time. The temperature of the 

E erturbation is derived from the rate of heating due to Thomson scattering of 
ackground radiation. Hutchins gives expressions for density, temperature, 
and electron concentration within the perturbation as functions of t^and-il.. 
Comparable calculations have been performed by Matsuda, Sato, and Takeda 
(1969) and by Lepp and Shull (1983). Results from these calculations may be 
used as initial conditions for studies of the later phases of cloud collapse. 

The collapse of gravitationally bound clouds in the post-recombination 
epoch has been studied by several investigators, all employing quite similar 
physical assumptions. In general the cloud has been taken as spherical and 
uniform with its evolution governed by (1) pressure-free Newtonian collapse; 

(2) H z formation by one or both of the reaction sets given above; (3) LTE 
excitation and optically thin emission from H z and H + ; and (4) energy balance 
employing spatially averaged temperature and density. Semiquantative 
descriptions of collapse based on comparison of free-fall, H 2 ~formation, and 
cooling times are given by Yoneyama (1972) and Silk (1977). Numerical 
integration of equations describing processes 1-4 has been performed by 
Saslaw and Zipoy (1967), Matsuda, et al. (1969), Hirasawa (1969), Hutchins 
(1976), Carlberg (1981), Palla et al (1983), and Lepp and Shull (1983). 

Some modifications to the general method have been the inclusion of a pressure 
gradient (Matsuda et al.) and the use of an escape probability (Carlberg) or 
Planck function (Palla et al.) for optically thick emission. 

These various studies have yielded qualatatively similar descriptions of 
cloud evolution, but considerably different predictions for the minimum Jeans 
mass reached during collapse. In clouds with initial temperature Id 2 - - 10 5 K 
and initial density 1-10^ cm , H z cooling keeps the temperature below 1500K 
until the density has increased to 10 IZ - cm -3 and the molecules are dissociated. 
During this evolution the Jeans mass falls by several orders of magnitude from 
10 r M to 1-200 M© , the minimum value depending on the detailed treatment of 
physical processes - especially H z cooling. After H t dissociation is com- 
pleted, the temperature rises nearly adiabatically and the Jeans mass increases 
until H is ionized at T^io 4 *, n-^ 10 ‘ <J .cm -3 . Absorption of ionization energy, 
together with free-free and bound-free emission, again slows the temperature 
increase, but as the cloud becomes optically thick the temperature resumes an 
almost adiabatic climb. The Jeans mass may have a second and smaller minimum 
during H ionization, but from the calculations performed to date there is no 
consensus as to whether the limiting mass for the first generation of stars 
is near that of the sun or 1-2 orders of magnitude larger. 

The present research was undertaken to extend and improve these studies 
of zero-metal clouds by treating in more detail the spatial variation of 
physical variables and the transfer of radiation in lines and continuum. With 
such a treatment it is possible to model the formation of a hydrostatic proto- 
stellar core and thus derive more directly a lower limit to stellar masses. 

In addition we investigate effects of fragmentaton on cloud evolution and 



protostellar mass. 


II. METHOD 
a) Hydrodynamics 

The hydrodynamic collapse of spherical gas clouds was simulated with 
a hydrodynamic computer code adapted "from the work of Bodenheimer (1968b). 
The differential equations of motion, mass continuity, and energy balance: 



are solved as implicit difference equations. The pressure P is that of an 
ideal gas, and the internal energy E includes terms for excitation of H* 
rotational and vibrational states, dissociation of H x and H~, and H ioni- 
zation. The form of the cooling function _/\_ is discussed in Section lie. 

Evolutionary model sequences have been constructed for clouds of mass 
10 s and 5 x 10 r M 0 , starting from a homogeneous, quiescent configuration of 
of density 4.8 x 10 _2?< gm/cm 3 and temperature 100 K. The chemical composition 
is taken to be 75 percent hydrogen and 25 percent' helium by mass, with all of 
the hydrogen initially in atomic form. The pressure at the cloud surface is 
assumed to remain at xts initial value. Spatial variation of the physical 
variables is represented with 100 spherical zones, and time steps are chosen 
so that changes in the quantity do not exceed 6 percent in a single step. 


b) Chemistry 

The numerical code includes formation of molecular hydrogen by the 
reactions (1) to (5) and also by three body combination involving helium: 

H + H + He — H l + He. ( 6 ) 

Molecular hydrogen is dissociated by the inverse reactions to (3) - (6), that 
is, by collision of H molecules with H z , H, and He. H atoms are ionized by 
collision with e~, H, H 2 , and He: 

H + e~ H + + e~ + e~ (7) 

H + H — > H - *" + e~ + H (8) 

H + Hz H + + e~ + Hz (9) 

H + He -=> H + + e~ + He, (10) 

and H+ recombines to H by- the inverse reactions to (7)- (10) (3-body recom- 
bination) and by radiative recombination: 

H + + e - — H + hy*. (11) 

Photoionizaion (the inverse of 11) is not considered because for most of 
their histories the clouds are optically thin to photons released in (11), 
and there is no external radiation field. Likewise, photodissociation is 



negligible. 

The reaction rates are identical to those used by Palla, et al., except 
that for reactions not considered by these authors (reactions 6 and its 
inverse; the inverses of 7 and 8; 9 and 10 and their inverses) rates were 
supplied by Prasad ( private communication ) . 


c) Cooling 

At densities less than 10 " lo gm/cm 3 , the radiative luminosity is 
primarily from rotational and vibrational transitions of collisionally 
excited H z . The tables given by Hollenbach and McKee (1979) are used for 
the local cooling rate _A_ as a function of the temperature and the H, H z , 
and He abundances. These calculations do not include excitation by He, but 
as the collisional cross-section of He is similar to that of H z (Green, et 
al. 1978), we use the combined abundance n(H 2 ) + n(He) in place of n(Hi). 

As the H t lines become optically thick, the cooling rate must be modified 
for trapping of photons within the cloud. We take 

A = ££n- A-.&E ... a (per unit volume), 

«l 0 * / 4j 

where n; is the population of the upper level, A;«* is the probability of a 
spontaneous transition from level j to i, &E • is the photon energy, and 
&•; is the photon escape probability. The sum is over the rotational 
transitions within the v=0 vibrational state and over vibrational transitions 
between v = 0, 1, and 2 for which the molecular parameters and escape proba- 
bilities are averaged over the rotational substates. For the optically 
thick transitions, the calculation of the nj is simplified by the fact that 
the density is large and the populations are in LTE. Also, calculation of 
the escape probability is simplified because the collapse velocity is much 
larger than the thermal velocity and ft>\\ may be determined from local physical 
parmeters by the method of Sobolev (1960). We take 
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where /U is the cosine of the direction angle 9 , V s is the fluid velocity 

in direction 9 , and h vs = A E,- • . 

y o 

At densities exceeding 10"'° gm/cm -3 , energy carried by the 0.75 ev 
photons produced in the reaction 

H + e~ — ^ H ~ + hv 

becomes the principal coolant. The numerical code treats emission and 
absorption of this radiation to determine the net local cooling . 

The emission coefficient is 

£ = h dn^- /dt = hv> (n H n e R H|e ) . 

The opacity for the H~ photons is mainly from free-free transitions of H~ 
(Chandrasekhar 1946; Ohmura and Ohmura 1960, 1961), H^ - (Somerville 1964), 
and He - (Somerville 1965). We use 
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If the cloud is optically thin to the H 
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radiation, the transfer equation 


in which I is the specific intensity, S =€/k is the source function, and 
dt^= “Kds is integrated along the outward and inward radial directions 
to obtain two intensities It and I~at each grid point. The mean intensity 
is obtained as a weighted mean 

J = u> I + + (1-<J) I~, 

where (P is the angular size of the small, high-density core in which the 
H~ photons are produced. The cooling rate is then obtained as 


_A = 4tr ~K (S - J) . 


If the cloud is optically 
diffusion approximation 

F = 


thick, the radiative flux is determined from the 
H7T AS 

~ T ar 


and the cooling rate from 
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III. RESULTS 


a) Evolution of a 5 x 10 5 " cloud 

Characteristic features of the evolution are illustrated in Figures 1 and 
2, which show the variation of temperature and Hj_ mass fraction with density 
at the cloud center. Initially, the cloud mass exceeds the Jeans massby a 
factor of 25 and the collapse Is nearly free-fall. During the evolution from 
points A to B the gas temperature increases almost adiabatically because there 
is insufficient H z for radiative cooling to balance compressional heating. 

(The quantity Y = d(ln T)/d(ln£> ) is actually slightly larger than 2/3 because 
of energy gained from the association of H. atoms.) However when the 
fraction has grown to 1.5 x 10~ H at point B, most of the compressional energy 
is radiated, and from B to C the temperature decreases. Eventually an 
equilibrium is reached with nearly equal timescales for heating, H z formation, 
and cooling, and between C and D J\_ it P d(l/ / o)/dt. At D the 3-body - 
reactions become highly effective, and there is increased heating from 
formation. Simultaneously, the cooling efficiency drops as the Ha. lines 
become optically thick. These effects cause the temperature to rise more 
steeply between D and E. At E the temperature has reached 2000 K and Hi 
begins to dissociate. The temperature increases less steeply from E to F as 
the compressional energy is used in dissociation. With the destruction of 
H z the gas can no longer radiate and from F to G there is again an adiabatic 
temperature rise. At point G hydrogen begins to ionize, and there is an 
abrupt increase in the rate at which free electrons react with hydrogen to 
form H~ . The escape of photons produced in this reaction keeps the 
temperature nearly constant from G to H as the density increases by 5 orders 
of magnitude, and because the increase in density favors Hj. formation over 
dissociation, the abundance attains a second and slighly higher maximum of 



X(H t ) =0.2. At point H the temperature again increases sharply as the gas 
becomes optically thick to the H ~ radiation. The temperature gradient is 
nearly adiabatic from this point until the calculation was ended at K, with 
only a slight moderation between I and J as H is ionized. 

Between points I and J the collapse decelerates near the cloud center 
and a small, nearly hydrostatic core Is formed. This condensation grows by 
accretion of infalling material and should eventually become a star (see, e.g. 
Stahler, et. al 1985). Its mass thus'represents the minimum mass of a star 
formed at the center of the cloud. Because there are no physical damping 
mechanisms, the core slightly overshoots its equilibrium configuration and 
begins a series of oscillations just before the calculation is terminated 
at K. At this point the core mass is 0.02 My and the accretion rate 
^TTr*- v just outside the core) is 0.6 M^/yr. 

If cooling by the photons produced with H ~ is neglected, the temperature 
rise would continue at point G In Figure 1, following the dotted curve. In 
this case the protostellar core has a mass of 0.5 M^, suggesting that only 
a massive star could be formed. 

The dashed curve in Figure 1 illustrates a collapse calculation by 
Palla, et al. (1983) that differs from ours in several important respects: 

(1) the density, temperature, and other physical parameters remain uniform; 

(2) the collapse is always free-fall; (3) the rate of optically-thick Hi 
cooling is estimated from a Planck function; and (4) cooling by H* photons 
is not considered. Although their numerical results resemble ours, the 
agreement beyond point D is only a fortuitous result of their more efficient 
Hi cooling compensating the absence of H~ cooling. Pala, et al. find that the 
Jeans mass decreases during the collapse to a minimum of 0.1 M G , but without 
spatial resolution they cannot identify the time at which a protostellar 
condensation forms. Our more detailed calculation provides a stricter lower 
limit for the stellar mass. 

The spatial variation of physical quantities within several cloud models 
is illustrated in Figure 3, which gives the variation of density, temperature, 
H 2 , abundance, and velocity with distance from the cloud center. The density 
is seen to be uniform near the cloud center and to decrease as r _z in the 
outer layers. This distribution is well known from previous studies of 
spherical collapse (e.g., Bodenheimer 1968a). The cloud collapses in free- 
fall until it is decelerated by a pressure wave moving inward from the 
surface. The uniform core is material not yet affected by the pressure 
gradient. The collapse velocity has a maximum between the edge of the core 
and the surface of the cloud. This maximum moves inward in radius and 
eventually steepens into a shock wave at the core boundary. The variation 
of temperature with radius is seen to resemble the variation of temperaure 
with density at the cloud center (Figure 1). This may be understood by 
noting that in the free-fall region the temperature is a function of density, 
and that in the more slowly evolving envelope the temperature at each point 
remains fixed at value related to the density of this point when it was at 
the core boundary. Likewise, the Hi abundance varies with radius in Figure 3 
much as it does with central density in Figure 2. 

b) Evolution of a 10^ M cloud 

The variation of central temperature and H z abundance with density in 
a 1C) 4 Mp cloud is shown in Figures 4 and 5. The mass of this cloud is near 
the minimum for instability to collapse. (In the initial configuration the 
mass is egual to half of the Jeans mass). The evolution starts as a series 
of adibatic oscillations between point A and A 1 until there is sufficient 
and radiative cooling for sustained collapse. Even during the collapse 
phase (until point D) the low-mass configuration is close to hydrostatic 
equilibrium. (Points in Figures 4 and 5 are labeled according to analagous 



points in Figures 1 and 2 . ) Again a balance is attained between the 
rates of compressional heating and radiative cooling, but these rates are 
smaller than in the case of free-fall collapse. The smaller cooling rate 
is reflected in the lower temperature between points B and D. However, the 
tracks of the low and high mass clouds converge when the 3-body reactions 
accelerate H formation and the cooling rate becomes more sensitive to 
density than temperature. 


c) Fragmentation sequences 

Two additional sequences of models were constructed to simulate the 
sepapration of a collapsing cloud into successive fragments of increasingly 
smaller mass. The times of fragmentation and fragment properties have been 
chosen in analogy with the numerical calculations of Bodenheimer (1978). At 
selected intervals the collapse is restarted from a uniform configuration at 
rest. It is assumed that fragmentation occurs when the ratio of the central 
density to the initial density corresponds to values given as Bodenheimer ' s 
Table 1, Squence ill; that the fragment mass is one eighth of the parent 
cloud mass; and that the fragment's density, temperature, and chemical 
composition have the values at the center of the parent cloud. 

The evolution of fragments of the 5 x 10^ M^cloud is illustrated in 
Figure 6 . The solid line shows the variation of central temperature with 
density in the parent cloud without fragmentation (Figure 1). The first 
fragment forms at point A, and its evolution follows the dotted line from A 
to B. Likewise, the evolutionary track of the second fragment begins at 
point B, and the third fragment at C. It is evident that the tracks of 
the first two fragments guickly approach and then closely follow the path 
of the parent cloud. This is because the compressional heating and 
radiative cooling are the same local functions of density and temperature 
in each case. The track of fragment C, however , _differs from that of the 
parent cloud in that the temperature begins a nearly adiabatic rise at a 
lower density. This is a result of the density homogenization at the time 
of fragmentation, which gives the fragment more high density material and 
a larger optical depth to H- photons, the principal coolant at this time. 

But as a more significance difference in the evolution with fragmentation, 
we note that fragment C is still in nearly free-fall collapse at the density 
where the parent cloud forms an equilibrium core. This lengthing of the time 
to reach equilibrium occurs because the size of the dense, collapsing core is 
increased with each fragmentation, and it takes longer for the core mass to 
shrink to the Jeans mass. We thus have not determined a minimum stellar mass 
in the case of fragmentation, but we note that as the Jeans mass is increasing 
at the point where the calculation was terminated, the initial mass of the 
protostellar core should be larger than in the case without fragmentation. 

Similar evolutionary tracks are obtained for fragments of the 10 M G 
cloud (Figure 7), although in this case the first fragment maintains a 
temperature higher than that of the parent cloud. The tracks differ because 
the fragment is collapsing more rapidly than the parent cloud, which is close 
to equilibrium, and the higher temperature increases cooling to balance the 
extra compressional heating. The tracks for the second and third fragments 
are almost identical to tracks for corresponding fragments in the 5 x 10 f M 0 
cloud . 
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EVOLUTION. OF ZERO-METAL CLOUDS 
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ABSTRACT. Dynamical, chemical and thermal evolution of zero-metal gas 
clouds has been modeled to study conditions of star formation in the 
early universe. Numerical results are shown for collapse of a spher- 
ical cloud of mass 5 x 10 ^ M®. Cooling by R 2 lines and by photons 
emitted in H + e- -*• H- + hv maintains collapse until formation of an 
equilibrium protostellar core with mass 0.03 M©. The cooling by 
photons produced with H- is essential for low-mass star formation. If 
the cloud is fragmented, the evolution of each piece resembles that of 
the parent cloud. 

1 . METHOD. 

Models for collapse of zero-metal clouds have been constructed with a 
numerical code based on the work of Bodenheimer 1968. The calculations 
include self-gravity, internal pressure, H 2 chemistry and radiative 
cooling with treatment of the spatial and time variation of all phys- 
ical quantities. Molecular hydrogen is formed by associative attach- 
ment of H with H- and by three-body combination. The gas is cooled by 
rotational and vibrational emission from H^, with a photon escape 
probability consistent with the supersonic collapse velocity, and by 
emission of 0.75 ev photons produced in H + e- -*■ H- + hv. The cooling 
associated with H- formation has not been considered in previous 
collapse studies. 

2. RESULTS AND CONCLUSIONS 

“•V. 

Figure 1 shows the evolution of central temperature with density in a 
5 x 10 5 M© cloud of 75 percent hydrogen and 25 percent helium which 
begins collapse from density 5 x 10 24 gm/cm^ and temperature 100 K. 

The temperature first increases sharply, reaching a maximum (point A) 
when sufficient H 2 is formed that the cooling time equals the compres- 
sional heating time. The temperature rise steepens again at R as the 
H 2 emission becomes optically thick and at C following H 2 dissociation. 
The onset of cooling by H- photons is abrupt at point D, and subsequent, 
evolution is nearly isothermal until the cloud is opaque to. this radi- 
ation at E. Without H- cooling the temperature rise would continue 
along the dotted curve. Point F marks hydrogen ionization and G the 


cessation of collapse in a small core. The core mass is 0.03 M© with 
the H- cooling and 0.05 M© without. Numerical accretion studies (e.g., 
Stahler et al. 1985) suggest, that the core mass will grow by a factor 
of at least 10 3 before the formation of a main sequence star. Thus 
only the models with H-_cooling could yield stars of low mass. 

The dashed curve in Figure 1 illustrates a calculation by Palla et 
al. for collapse and cooling of a pure H cloud. Although- their numer- 
ical results resemble ours, the agreement is likely fortuitous as the 
physical basis of cooling differs in the two calculations. The Palla 
et al. model does not include H- cooling, but it has more efficient 
cooling in the optically thick Hj lines. 

If our cloud model is assumed to fragment during collapse, the evo- 
lutionary tracks of the fragments quickly converge to that of the parent 
cloud. Fragmentation may thus provide multiple centers for collapse, 
but it should not affect the stellar mass function. 



Figure 1. Temperature and density evolution in zero-metal clouds. 
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